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Computer programming before 1970 


= Elegancy: not a main issue 
= No general methodology 
= No structure 


= Main computer languages: machine code, 
assembler, FORTRAN, COBOL 


= Main goal: Fit the program into memory 


A computer program 





Structured programming appears 

= E. W. Dijkstra, Go to statements considered 
harmful (1968) 

= E. W. Dijkstra, Notes on structured 
programming (1972) 

= D. E. Knuth, Structured Programming with 
Go-To Statements (1974) 


Structured programming spreads 


= Kosaroju (1974), D-charts 
= McCabe (1976), representation using graphs 


= Appliying structured programming to the 
formulation of algorithms 
— Henderson and Snowdon (1972) 
— Knuth and Szwarcfiter (1974) 


Structured programming takes over 


With the time became a general technique for the 
formulation of algorithms 


A natural question 


= How to tell if a given algorithm is structured ? 
= Apparently not yet considered. 


m= How to tell if two structured algorithms have a similar 
structure (syntax) ? 


= Apparently not yet considered. 


Our proposal 


= Introduce a class of graphs which corresponds 
to structured algorithms, named Dijkstra graphs 


= Answer the question of recognizing structured 
algorithms (i.e. Dijkstra graphs). 


= Answer the question whether two structured 
algorithms are similar (i.e. solving the 
isomorphism problem for Dijkstra graphs). 


The proposed algorithms for recognition and iso- 
morphism of Dijkstra graphs have both linear time 
complexity. 
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Basis of Structured Programming 


E. W. Dijkstra, Notes on structured programming, in 
Structured Programming, Academic Press, 1-82, 1982. 


Three types of statements: 


Sequence: 


: 


a P 


Basis of Structured Programming 


Selection 


SELECTION 


IF C THEN P 


IF C THEN P, ELSE P, 





Basis of Structured Programming 


Selection: 


CASE C OF P,, ... Px 


Basis of Structured Programming 


Iteration 


T 
27 i) WHILE C DO P 
F 
: REPEAT P UNTIL C 


Our Problem 


Given a computer program (flowcharts) Is it 
structured ? 


Approach: graph theory 


Flow graphs and reducible graphs 


= A flow graph is a directed graph G containing 
One source s(G) 
= A reducible graph is a flow graph G such that 
every cycle C’ contains a unique entry point, for 
paths starting in s(G) 
s(G) 


Flow graphs and reducible graphs 


m= Every algorithm (structured or not) corresponds 
to a flow graph. 


m= Every structured algorithm corresponds to a 
reducible graphs. 


Reducible graphs 


= Characterizations: Hecht and Ullman (1970, 
1974), Tarjan (1974) 


= Efficient Recognition: Tarjan (1974) 


= Useful characterization: A graph G is reducible 
if and only if any DFS of G, starting from s(G) 
determines the same set of cycle edges. 


s(G) 


cycle 
edge 


Statement graphs 


A statement graph Is a directed graph H, satisfying 

= Each vertex is labelled as regular (R) or 
expansive (X) 

= [7 belongs to one of the following classes: 

trivial graph 

sequence graph 

if graph 

if-then-else graph 

— p-case graph, p> 3 

while graph 

— repeat graph 


Statement graphs 


(a) (b) (Cc) 


Note: All single source-sink 


xX 


Statement graphs 
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(c) (f) (g) 


Note: All single source-sink 


Closed subgraphs 


Let H, single source-sink subgraph of G 
H closed when 

1.vEV(A)\ s(f1) > N-|v| CV(A) 
2.ve€V(A)\t(A) > N*|v] CV(A) 

3. us(H) is a cycle edge => v € N*(s(#)) 


s(H) 


t(H) 


Prime subgraphs 


Let H be an induced subgraph of G 
His prime when #4 is 


1. isomorphic to some non-trivial statement graph 


2. closed 


NOTE: The while and repeat graphs might be iso- 
morphic as statement graphs, but distinguishable as 
prime subgraphs, even with no labels. 


Expansion operation 


Let G, H graphs, V(G) NV (H) = 9, H single 
source-sink, v € V(G). 
The expansion of v into H consists of replacing uv by 
H, in G, such that 

= Na(s(H)) = Nav) 

= NG (t(H)) = Nev) 

= Remaining adjacencies unchanged 





Contraction operation 


G graph, H prime subgraph of G 
The contraction of H into a single vertex: 


= Identifies the vertices if H into the source s(/) 
of H 


= Remove possible parallel edges or loops. 





Dijkstra graph 


A A Dijkstra graph (DG) is a graph with vertices 
labelled as X or R, recursively defined as: 


1. A trivial statement graph is a DG 


2. Any graph obtained from a DG, by expanding 
some X-vertex of it into a statement graph is 
also a DG. Furthermore, after expanding an 
X -labelled vertex into H, the vertex s(/) is 
labelled as R. 


Kxample 





Example: Binary search 


void main() 


{ 
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Kt 
ro 
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ro 
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int c, first, last, middle, n, search, array[100]; 
scanf ("%d",&n); //number of elements 
eG = 05 
while (c < n){ 
scanf ("Sd", &array[c]); 
Cres 


} 


scanf ("$d", &Search) ; 





//value to find 


first = 0; last =n - 1; middle = (firsttlast) /2; 
while( array[middle] != search && first <= last ){ 
1f ( array[middle] < search ) { 
first = middle + 1; 
middle = (first + last) /2;} 
else{ 
last = middle - 1; 
middle = (first + last) /2;} 
} 
if (item == a[mid]) { 





printf ("Binary search successfull!!\n") } 
else { 


printf("Search failed! \n") } 


Constructing Dijkstra graphs 


Theorem 1 A graph G is a DG if and only if there is 

a sequence of graphs Go,...,G, each of them 

having the vertices labelled as X or R, such that 

1. Go Is the trivial graph, with the vertex labelled X 

2.G,=G 

3. G; is obtained from G;_1, 1 > 1, by expanding 
some X-vertex of it into a non-trivial statement 
graph, whose source s(H) receives label R. 


Still: How to recognize ? 


Basic properties 


Lemma 1 Let G be a DG. Then 

1. G contains some prime subgraph 
2. Gis single a source-sink graph 
3. Gis a reducible flow graph 


Independence 
H(G') = set of non-trivial prime subgraphs of G, 
Let H, H’ ¢ H. H, H’ are independent when: 


V(A)OV(A') = 9, or 
V(H)OV(#’) = {v}, where 
v= s(H) =t(HA’), orv = s(A’) = t(#) 


sf) 


S(A’ 


1(H") 
1(H") 


Independent primes 


Lemma 2 Let H, H' be two distinct prime 
subgraphs of some graph G. Then H, H’ are 
independent. 


Images 


G graph, 

H(G') = set of non-trivial prime subgraphs of G, 
H €H(G) 

G | H = graph obtained from G by contracting H 


Image of a vertex: 
For v € V(G), the image of v inG | H = 


fante)= 4° ,ifu gé V(A) 


s(H)_ , otherwise. 


Images 


Image of a subset of vertices: 
For V’ C V(G), the (subset) image of V’ in G | H = 
Teja (V") = Wvew'Leya(v) 


Image of subgraph: 
For H’ C G, the (subgraph) image of H’ inG | H = 


Ig,y(’) = subgraph induced in G | H by the sub- 
set of vertices Ie, 7(V(7’)) 


Preservation of Primes 


G, arbitrary graph 
H,H' €H(G), HF H". 


Lemma 3 [cq 47(H') © H(G | A) 


Commutative law 


Lemma 4 H, H’ € H(G) 


(G1 A)l Uqua(#")) = (G1 A") 1 daa (7) 


Contractile sequences 
A sequence of graphs Go,..., G; iS a contractile 
sequence for a graph G, when 


1.G=Go, 


2. Gis, = G; | Aj, for some A; € H(G;), 7 < k. Call 
H,, the contracting prime of G;. 


Go,..., Gz iS maximal when H(G;,) = @. 


In particular, if G; is the trivial graph then Go,..., Gx 
IS maximal. 


Iterated images 


Go,...,G,, contractile sequence of G 
H,, contracting prime of G; 
(1.e. Giy41 = G; al H;) 


For H; C Gj, the iterated image of H; in G,, q => 7 is 


A, ifg=7 
Ig,Ue,.,(H;)), otherwise. 


Ig,(H;) = 


Main characterization 


Theorem 2 Let G be an arbitrary graph, with 
Go,...,G, and Go,...,G,, two contractile 
sequences of G. Then G; and G‘,, are isomorphic. 
Furthermore, k = k’. 


Notes: 
1) Leads to a greedy recognition algorithm. 


2) The input graph G has no labels. 
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iterated image of H, 


Gy41 = Gy ts and Gray = Gy, i; 


Go— Gi, Gi = Gy n~G: =] G,Ga7 CG a= 
Hg = Hod = Aigesigdl ge 1, 


Recognition 


Corollary 1 LetG be an arbitrary graph, and 
Go,...,G, a contractile sequence of it. ThenG Is a 
DG if and only if G;, Is a trivial graph. 


Bottom-up contractile sequences 


= G, reducible graph 
= Go,...,G,, contractile sequence C of G 
= /; the contracting prime of G;,0<i<k 


= C IS a bottom-up (contractile) sequence of G when 
each contracting prime H; satisfies: 
s(H;) is not a descendant of s(#), for any prime 


Generating primes by contractions 


Contracting a prime H of G may generate a (new) 
prime H’ of G| H. 


H' 


Generator and generated 


Relation between generator and generated primes: 


Lemma 5 G, reducible graph 
HeEH(G), He H(GIlA)\H(G),v=s(A). 
Then v is a descendant of s(H') inG | H. 


Recognition Algorithm 


= G, reducible graph 

= lteratively, find a lowest vertex v of G, s.t. 

= v Is the source of a prime subgraph H of G 
m= then contract H. 

= otop when no primes exist any more. 


The initial graph is DG iff the final is trivial 


The algorithm 


G, arbitrary flow graph (no labels) 
Ec, set of cycle edges of a DFS of G, starting at 
s(G) 
U,-.+,Un, topological sorting of G — Ec 
=7 
while 7 > 1 do 
if G is the trivial graph 
then return YES, stop 
if vu; is the source of a prime subgraph A of G 
thenG :=G fA 
1:=1—1 


return NO, stop 


Avoiding to recognize reducible graphs 


Lemma 6 LetG be an arbitrary flow graph. If the 
algorithm recognizes G as a Dijkstra graph then 
necessarily G is a reducible graph. 


Complexity 


= lopological sorting: O(m) 


= Deciding if vertex uv; is the source of a prime: 
O(|N*(u;)]) 


= Ihere can be O(n) primes 
= Contracting a prime: O(|.N*(v;)|) 


m Each edge is traversed only a constant number 
of times, overall 


= Complexity: O(m) 


Bounding the size of a DG 


Lemma 7 LetG be a DG graph. 
Then m < 2n — 2. 


Consequence: 
Complexity of the recognition algorithm: 


O(n) 





void main() 


{ 
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int c, first, last, middle, n, search, array[100]; 
scanf ("Sd",&n); //number of elements 
c= 0; 
while (c < n){ 
scanf("Sd",éarray[c]); 
ctt; 
} 
scanf ("%d", &search); //value to find 
first = 0; last =n - 1; middle = (firsttlast)/2; 
while( array[middle] != search && first <= last ){ 
if ( array[middle] < search ){ 
first = middle + 1; 


middle = (first + last) /2;} 
else{ 
last = middle - 1; 
middle = (first + last) /2;} 
} 
if (item == a[mid]) { 


printf ("Binary search successfull!!\n") } 
else { 
printf ("Search failed! \n") } 


void main() 


{ 
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22 } 


int c, first, last, middle, n, search, array[100]; 
scanf ("%d",&n); //number of elements 
c= 0; 
while (c < n){ 
scanf ("Sd", éarray[c]); 
Cots 
} 
scanf ("3d",&search); //value to find 
first = 0; last =n - 1; middle = (firsttlast) /2; 
while( array[middle] != search && first <= last ){ 
if ( array[middle] < search ) { 
first = middle + 1; 


middle = (first + last) /2;} 
else{ 
last = middle - 1; 
middle = (first + last) /2;} 
} 
if (item == al[mid]) { 


printf ("Binary search successfull!!\n") } 
else { 
printf("Search failed! \n")} 





void main () 


{ 


DAYANA PF WNE 


be wo 
oO 


11 
12 
13 
14 
15 
16 
17 
18 
18 
20 
2d 
22. 3 


int c, first, last, middle, n, search, array[100]; 
scanf("%d",é&n); //number of elements 
c= 0; 
while (c < n){ 
scanf("%Sd",éarray[c]); 
ctt; 
} 
scanf("%d",&search); //value to find 
first = 0; last =n - 1; middle = (firsttlast)/2; 
while( array[middle] != search && first <= last ){ 
if ( array[middle] < search ){ 
first = middle + 1; 


middle = (first + last) /2;} 
else{ 

last = middle - 1; 

middle = (first + last) /2;} 


void main() 


{ 
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22 } 


int c, first, last, middle, n, search, array[100]; 
scanf ("$d",&n); //number of elements 
c= 0; 
while (c < n){ 
scanf ("Sd", &array[c]); 
Orr; 
} 
scanf("%d",&search); //value to find 
first = 0; last =n - 1; middle = (firsttlast) /2; 
while( array[middle] != search && first <= last ){ 
if ( array[middle] < search ) { 
first = middle + 1; 


middle = (first. + Last) /27} 
else{ 

last = middle - 1; 

middle = (first + last)/2;} 
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void main () 


{ 
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22 } 


int c, first, last, middle, n, search, array[100]; 
scanf("Sd", én); //number of elements 
c= 0; 
while (c < n){ 
scanf ("sd", é&array[c]); 
ctt+; 
} 
scanf ("$d",ésearch); //value to find 
first = 0; last =n - 1; middle = (firstt+last)/2; 
while( array[middle] != search && first <= last ){ 


void main () 


{ 
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22 } 


int c, first, last, middle, n, search, array[100]; 
scanf ("$d",&n); //number of elements 
c= 0; 
while (c < n){ 
scanf ("%Sd", éarray[c]); 
Gk 
} 
scanf ("%d", &search); //value to find 
first = 0; last =n - 1; middle = (firsttlast)/2; 
while( array[middle] != search && first <= last ){ 
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void main () 


{ 
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22} 


int c, first, last, middle, n, search, array[100]; 
scanf ("%d", én); //number of elements 
c= 0; 
while (c < n){ 
scanf ("%Sd", éarray[c]); 
Cry. 
} 
scanf ("%d",&search); //value to find 
first = 0; last =n - 1; middle = (firsttlast) /2; 


void main () 


{ 
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22 } 


int c, first, last, middle, n, search, array[100]; 
scanf ("$d",&n); //number of elements 
c= 0; 
while (c < n){ 
scanf ("%Sd",éarray[c]); 
Grsak 
} 
scanf ("%d",&search); //value to find 
first = 0; last =n - 1; middle = (firsttlast)/2; 
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void main () 


{ 


int c, first, last, middle, n, search, 
scanf ("%d", én); //number of elements 
c= 0; 
while (c < n){ 

scanf ("%sd", &array[c]); 

ctt; 


void main () 


{ 


int c, first, last, middle, n, search, 
scanf("Sd",&n); //number of elements 
c= 0; 
while (c < n){ 

scanf ("%d",é&array[c]); 

Cr+: 


array[100]; 


array[100]; 


22 


ze 


void main () 


scanf ("Sd", &n); 
c= 0; 


void main () 


scanf ("%d", &n) ; 
c = 0; 


int c, first, last, middle, n, search, 


//number of elements 


int c, first, last, middle, n, search, 


//number of elements 


array[100]; 


array[100]; 


OorArnuA WHE 


void main() 


{ 


Isomorphism of DG graphs 


General idea: 


= Define a code CG) for a DG graph G, 
such that C(G) uniquely identifies G. 


a G, = G, iff C(G,) = C(G2) 
= ((G) is a string of integers C {1,..., A(G) + 4} 


Coding a DG graph 


= Assign an integer, named type(H), for each statement 
graph H, 
= acode C'(v) for each vertex v € V(G), 


= acode C(#H) for each prime subgraph H of a bottom-up 
contractile sequence of G. 


= Ihe code C(G) of the graph is defined as 
C(G) := C(s(G)). 

= For a subset V’ C V(G), the code CV’) of V’ is the set of 
strings C(V") = {C(u;) |v; € VF. 

a Write lex(C(V"’)) = C(w,)|I...||C(w,), 
whenever V’ = {v,,...,v,} and C(v;) is lexicographically 
not greater than C'(v;.1). 


type(H); C(v), v=s(H); CCH) 


a 


4\|O(N+(v) NN~ (W)C @) \ N= (2) 
SEWN FWIIICWHP@) \ fo) 
__itthen-else | 6 | 6lldew(CIW (NIC?) 

p+ Alitex(C(N*(»)) ICN?) 





Setting the codes 


= |he types of statement graphs: table 
= Encoding C(wv), v € V(G): initially set to 1. 


= Subsequently, if v becomes the source of a 
prime graph H, assign C(v) := C(v)||C(A) 


= Encoding C(#): table 


= ((/) is written in terms of type( 1) and 
C(w)|w € V(#), and so iteratively. 


Isomorphism algorithm 


G, DG; Eco, set of cycle edges 


Find a topological sorting v1,...,Un of G— Eo 
fort =n,n—1,...,1 do 
C'(v;) =I 


if v; is the source of a prime subgraph H then 
2||C(N* (v;)), if H is a sequence graph: 
3||C(NT (v;) \ NT2(v;))||C(NT?7(v;)), if A is an if-then graph; 
AIIC(NT (vi) ON (v:)) ICON (vi) \ N7 (i), 

if AT is a while graph, 
5|IC(N* (vi) ||O(N** (va) \ {vi}, 

if H is a repeat graph; 
6||lex(C(N*(v4)))||C(N*? (u:)), 

if AZ is an if-then-else graph. 
p + Al|lex(C(NF (vi)))|IC(NT* (vi), 

if H is a p-case graph. 


Clu;) = C(vi) || 


C(G) := C(v1) 


Kxample 





G, Dijkstra graph 


Kxample 





C'(v19) _— 12||C'(vj44) = 121 


Kxample 








C(v9) = 16||lea(C(v11), C(vi2))||C(vig) = 16111 


Kxample 








C(v4) = 12||C(vg) = 1216111 


Kxample 





Kxample 





C(ug) = 12||C(v7) = 121 


Kxample 








C'(v3) = 13||C(vg)||C(vg) = 131211 


Kxample 








Cv) = 14||C(vs)||C(vs) = 1412161111 


Example 








Cv) = 16||lex(C(v2), C(vs))||C(ro) 
C(v,) = 161312111412161111121 


Example 


C(v,) = 161312111412161111121 


Correctness 


Theorem 3 Let G, G’ de Dijkstra graphs, and 
C(G), C(G’) their encodings, respectively. Then 
G,G' are isomorphic if and only if C(G) = C(G’). 


Consequences 


Corollary 2 Let G be a DG and C(G) its encoding. 


1. There is a one-to-one correspondence between the 1's of 
C'(G) and the vertices of G. 


2. The encoding C(G) of G is unique and is a representation 
of G. 


Exhibiting the isomorphism function: 

Corollary 3 LetG,G’ be DGs and C(G), C(G’) their 
corresponding encodings, satisfying C(G) = C(G’). Then an 
isomorphism function f between G and G’ can be determined 
as follows. Letv € V(G) andv' € V(G") correspond to 1'°s at 
identical relative positions in C(G) and C’(G"), respectively. 
Define f(v) :=v'. 


Complexity 
Lemma 8 Let G be a Dijkstra graph, and C(G) its 
encoding; Then |C(G)| < 2n — 1. 


Theorem 4 For DGs, the isomorphism algorithm 
terminates within O(n) time. 


Generalization 


More useful statement graphs. 
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BREAK-WHILE BREAK-REPEAT 
GRAPH GRAPH 


Generalization 





CONTINUE-WHILE CONTINUE-REPEAT 
GRAPH GRAPH 


Applications 


Recognition algorithm: 
Graph watermarks 


lsomorphism algorithm: 
Code similarity analysis 


THE END 
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